Polarization effects in diffraction of light on a planar chiral structure 
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We analyzed polarization changes of light diffracted on planar chiral array from the standpoint of 
the Lorentz reciprocity lemma and found the bi-orthogonality in the polarization eigenstates for the 
waves diffracting though the grating in the opposite direction. Both reciprocal and non-reciprocal 
component in the polarization azimuth rotation of the diffracted light were identified. The structural 
chirality of the array arrangement and the chirality of individual elements of the array give rise to 
polarization effects. 
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Recently we reported that planar (2D) chiral struc- 
tures affect the polarization state of light in an enan- 
tiomeric fashion, similarly to three-dimensional chiral 
media lj. However polarization phenomena on diffrac- 
tion from planar chiral structures have never been stud- 
ied theoretically before, leaving the fundamental prop- 
erties of 2D chirality not fully understood. Here we re- 
port on the results of a theoretical investigation of polar- 
ization changes for light diffracted on regular arrays of 
planar chiral metallic structures from the standpoint of 
the Lorentz reciprocity theorem. By analyzing propaga- 
tion of light in two opposite directions we have identified 
a strong component in polarization effect on diffraction 
that can be induced either by the chirality of the indi- 
vidual elements of the array, or by arranging non-chiral 
elements of an array in a chiral fashion. 

Let us consider a planar square periodic array of metal- 
lic elements of thickness t with equal pitch d along the 
axes x and y placed between planes z = and z = — t 
(see Fig. ^and Fig. EI> • If a plane electromagnetic wave 
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of unit amplitude and polarization vector Aj is incident 
on the array from the region corresponding to z > 0, the 
transmitted field may be written as a summation over all 
diffracted waves, numbered by integer indices q and p 
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where & qp and h qp are amplitudes and wave- vectors par- 
tial of diffracted waves and 



k gP = g + h gP - e z \/k 2 - |g + h gp | 2 , 



h qp = 2Tr(qe x +pe y )/d. 



Here e x , ej, and e z are unit vectors along the axes x, y 
and z, g is the component of transverse to the axis 
z and k — |kj|. Let us now consider a "reversed" wave 
with polarization vector A r approaching the array from 



the opposite side of the structure (z < —t) along the 
direction of one of the partial diffracted waves of the " di- 
rect" scenario, with indices s and I and the wave vector 

kr = k s £ , 
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In the region z > this wave will produce diffracted 
waves with amplitudes h qp . This corresponds to the "re- 
versed" scenario of diffraction. The Lorentz reciprocity 
lemma 2] applied to the field superposition in the volume 
bounded by surface S that consists of planes x = ±d/2, 
y = id/2, z — z\ > 0, and z — Z2 < — t may be written 
in the following form: 



{[Ei x H r ] - [E r x isAdtr = 



(4) 



where Ej, H; and E r , H r are electric and magnetic fields 
created by waves incident from opposite directions. By 
using corresponding field expressions it may be shown 
from formula (4) that 



V^ 2 - |g| 2 (Ai • b sl ) = Vk 2 - |g + h si | 2 (A r • a si ). (5) 

The equality (0 constitutes the universal relation be- 
tween the amplitudes of partial waves in the " direct" and 
"reversed" diffraction scenarios. Scattering processes are 
often described in terms of 2 x 2 transformation matrices, 
relating Cartesian components of electric fields in coor- 
dinate frames of incident and scattered waves. For the 
"direct" (D) and "reversed" (R) scenarios these matri- 
ces for the incident and partial diffracted waves can be 
introduced as follows: & s i = DAi, h s i = RA r . It may be 
shown from equation (5) that these matrices are linearly 
related and mutually transposed: 
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index and c = 
\Jk 2 — |g + h s i\ 2 /y^k 2 — |g| 2 . For the purpose of 
analyzing the polarization eigenstates of the diffraction 
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FIG. 1: Coordinate systems and waves in the "direct" (a) and 
"reversed" (b) diffraction scenarios. 



process it is instructive to present both scattering ma- 
trices in the coordinate frame of the direct scenario (see 
Fig. l,a) where the operator of the reversed scattering 
process acts on the complex-conjugated field amplitudes. 
Here the polarization eigenstates are simply two linearly 
independent eigenvectors of matrices R* and D. The 
relation between them may be derived from equation 
©, which, when converted to the coordinate frame 
of the direct scenario, gives R* lm — cD* mn . It follows 
from the theory of matrix operators that eigenvectors 
of the Hermitian-conjugated matrices with elements 
D nm and D* nn , and therefore of matrices R* and D, 
are biorthogonal or, in terms of polarization eigenstates, 
are represented by antipode points on the Poincare 
Sphere, as shown on Fig. 3. In general, the point 
representing the first eigenstate in the direct scenario 
Id is an antipode to one of the points, which represents 
the eigenstates of the reversed scenario (this point is 
designed as 2r in Fig. 3). However eigenstate Id does 
not necessary coincide with eigensate lr of the reversed 
scenario, nor eigenstate 2d coincide with 2r. Therefore, 
the polarization eigenstates in the direct and reversed 
scenarios presented in the coordinate frame of the direct 
scenario could be different. Such a situation takes place 
if the complex matrix D is an asymmetric or even a 
non-diagonal matrix. 

We found that scattering matrices of non-zero order 
diffraction on periodic plano-chiral arrays, where chiral- 
ity is due to either structural chirality, or chirality of 





FIG. 2: Structural elements of the arrays: planar straight 
cross tilted against the array greed on the tilt angle if>, and 
chiral right-handed gammadion with bending angle <p. 
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FIG. 3: Schematic representation of diffraction on the 
Poincare sphere. For a chiral grating polarization eigenstates 
in the "direct" scenario (Id and 2d) and "reversed" scenario 
(lr and 2r) are elliptical. Polarization eigenstates (under- 
lined) for non-chiral gratings are mutually perpendicular lin- 
ear polarization. 



individual elements of the array, are either asymmetric 
or non-diagonal. Below we will illustrate these properties 
by numerical modelling of the diffraction process for var- 
ious planar chiral gratings. We calculated the fields and 
polarization characteristics of light diffracted on gratings 
numerically using the method described in reference |3j. 
It is based on a vector integral equation for the surface 
current induced by the light wave on the array particles. 
The equation is derived with boundary conditions for 
ideal metallic structures that assume a zero value for the 
tangential component of the electric filed on the metal. 
The integral equation is then reduced to an algebraic 
equation set by use of the Galerkin technique. 

In our modelling we concentrated on planar chiral ar- 
rays of 442 symmetry wallpaper group and calculated 
the polarization eigenstates of the diffraction process and 
polarization changes occurring in the diffracted wave for 
different incident polarizations. 

We studied diffraction for two different incident angles 
f3, at Pi = and at fa = arcsin(7r/fcri). In the first case, 
the diffracted wave (q = 1, p = 0) propagates at an- 
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TABLE I: Polarization Eigenstates (PES) for various diffraction processes presented in the "direct" scenario coordinate frame 
(all angles are measured in degrees, subscripts d and r denote "direct and "reversed" scenarios) 



Structure 


"Direct" Scenario 


"Reversed" Scenario 


Type of diffraction 


1-st PES 
(degrees) 


2-nd PES 
(degrees) 


1-st PES 
(degrees) 


2-nd PES 
(degrees) 


Straight crosses 
ip = 0; f5 = 


On = o.oo 
Vid = o.oo 


6»2d = 90.00 
r?2d = 0.00 


Olr = 0.00 

rn r = 0.00 


2r = 90.00 

V2r = 0.00 


No chiral effect 

Old = Olr, Vid — Vlr 


Straight crosses 
ip = +15; = 


Qi d = -14.4 
r, ld = -0.05 


6>2d = 77.1 
r?2d = 0.08 


Olr = -12.9 
rj lr = -0.08 


2r = 75.6 
?72r- = 0.05 


Chiral effect is present 

Old — Olr = —1.5, #2d — ^2j- = 1.5 


Straight crosses 
ip = -15; = 


Sid = 14.4 
Jftd = 0.05 


e 2 d = -77.i 

r/ 2d = -0.08 


6>i r = 12.9 
Vlr = 0.08 


Olr = -75.6 

T]2r = -0.05 


Chiral effect is present 

Old — Olr = 1.5, #2d — $2r = —1.5 


Right gammadions 
cp = 120; = 


fid — 0.1 

r?id = 4.38 


(72d — — ^0.0 

r/2d = -0.02 


"It- — Do.O 

r?i r = 0.12 


U2r — — OO.y 
T]2r = -4.42 


Chiral effect is present 
Old - 0i r = -57.4, 6»2d - Oar = 57.3 


Left gammadions 
if = 120; = 


fid = -6.1 
Vid = -4.38 


6»2d = 26.6 
V2d = 0.02 


0i r = -63.5 
rur = -0.12 


2r = 83.9 

7]2r = 4.42 


Chiral effect is present 
On - fir = 57.4, 6»2d - 02r = -57.3 


Right gammadions 
ip = 120; = 2 


Bid = 9.9 
?7id = 5.2 


y 2d = -80.0 
ri2d = 5.2 


Olr = 9.9 
?7lr = —5.2 


6>2r = -80.0 
r/2r = —5.2 


Chiral effect is present 

Old = Olr, Vid = ~Vlr 


Left gammadions 
ip = 120; = 2 


9 ld = -9.9 
»?id = -5.2 


6»2d = 80.0 
r\id = -5.2 


Olr = -9.9 
r,i r = 5.2 


6»2r = 80.0 
r; 2r = 5.2 


Chiral effect is present 

Old = Olr, Vid — — Vlr 



gle £ = £1 = arcsin(27r/fcc?) to the array. In the second 
case the diffracted wave (q = — 1, p = 0) the same angle 
makes to the array as the incident wave £ = £2 = — $2 
(for definitions of angles see Fig.l). The wave's polariza- 
tion azimuth 9 and degree of ellipticity 77 were calculated 
from the Cartesian field amplitudes using the standard 
definitions: tan2# = S2/S1, sin2?7 = s^/so, where Sj are 
the Stokes parameters. The results of our analysis for 
d = 4 /im, A = 2ir/k = 1520 nm, ft = 0, & = 22.3°, 
02 = 11.0°, £2 = —11.0° are summarized in Tableland 
Fig.0] We considered an array without a substrate. The 
width of the metal strips was equal to 0.05 /im. 

In optics polarization elements are often classified as 
reciprocal or non-reciprocal depending on whether their 
effect on the polarization state of the transmitted light is 
the same or different for light propagating in the oppo- 
site directions. This understanding of optical reciprocity 
which we will use below is somewhat different from the 
general, more tolerant definition of reciprocity based on 
the Lorentz lemma. For the purpose of comparison of 
the polarization transformations for opposite directions 
of light propagation, in the table and figures the polar- 
ization parameter of the waves are converted into the 
coordinate frame of the direct scenario. In such a pre- 
sentation if the values of polarization azimuth rotation 
in the direct and reversed scenarios are the same, the ro- 
tation is truly non-reciprocal. On contrary, a difference 
between the values of polarization azimuth rotation in 
the direct and reversed scenarios would represent a re- 
ciprocal component of the polarization change, that is 
analogous to the optical activity effect in a chiral liquid. 

The calculations revealed that: 



i) For all diffraction processes involving twisted or non- 
twisted arrays, equalities (6) are held to within the nu- 
merical accuracy of the method. They are thus compat- 
ible with the Lorentz lemma. 

ii) For the arrays of straight crosses polarization az- 
imuth rotation in opposite directions have opposite signs 
due to the difference in the efficiency of diffraction 
for perpendicular polarization components (line C on 
Fig. 4(B)). This non-reciprocity of polarization azimuth 
rotation is analogous to the polarization rotation non- 
reciprocity in dichroic media due to anizotropic dissipa- 
tion. 

iii) No polarization rotation is seen in the non- 
diffracted part of the beam at the normal incidence. It's 
polarization eigenstates are the same in both directions 
and for any type of array. 

iv) From Table I, one can see that for non-zero order 
asymmetrical diffraction (| £ |^| (3 |) when individual 
structural elements of the array are twisted rosettes po- 
larization azimuths of eigenstates for the "direct" and 
"reversed" scenarios are resolutely different. The dif- 
ference between the polarization azimuths of the eigen- 
satates depends on the rosette curvature angle (f and 
reaches a maximum of about 57° at ip = 120°. The 
difference in the eigenstates vanishes at rosette bending 
angle ip = 95°. 

v) When individual structural elements of the array 
are twisted rosettes, polarization azimuth rotation on 
diffraction has both reciprocal and nonreciprocal compo- 
nents. The non-reciprocal component of the polarization 
azimuth rotation is due to a difference in the efficiency 
of diffraction for perpendicular polarization components. 
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FIG. 4: Polarization azimuth rotation A = 9 — a on diffrac- 
tion from chiral arrays as a function of incident polarization 
azimuth, straight line - "direct" scenarios, dashed line - "re- 
versed" scenarios (all results are presented in the direct sce- 
nario coordinate frame). The influence of chirality is mani- 
fested as a split between corresponding solid and dashed lines: 
A) Array of straight crosses. The insert shows the chiral dif- 
ference in polarization azimuth for the two polarization eigen- 
states as a function of the tilt angle of the crosses; B) Array 
of left (L) and right (R) gammadions. The inset shows the 
chiral difference in polarization azimuth for the two polariza- 
tion eigenstates as a function of gammadion bending angle. 



The corresponding oscillating dependence of the nonre- 
ciprocal rotation on the incident angle is shifted in re- 
spect of line C corresponding to straight crosses. It is 
shifted along the incident polarization azimuth axis to- 



wards left for left rosettes and towards right for right 
rosettes. The split between corresponding solid and 
dashed lines in Fig. 4 indicates the reciprocal component 
of the polarization azimuth rotation analogous to optical 
activity. 

vi) Nonreciprocity of polarization rotation in the 
diffraction process is evident when a diffracted light wave 
is reflected straight back towards the twisted planar 
structure by a mirror, and then diffracts again. The po- 
larization state of the returning light after the second 
diffraction is different from that of the incident light, 
even if the incident light was an eigenstate in the for- 
ward direction. For an array of rosettes with ip = 120°, 
the two incident eigenstates and corresponding returning 
polarizations have azimuths different by 27° and 21°. 

Therefore, polarization effects on diffraction from pla- 
nar chiral grating can be induced by either structural 
chirality or the chirality of individual elements of the ar- 
ray. However, no polarization rotation compatible with 
the Lorentz lemma is possible for a wave transmitting 
through or reflected from a planar chiral structure at nor- 
mal incidence, as the scattering matrices are diagonal in 
this case. 

Finally, we shall note that our analysis is underpinned 
by the Lorentz lemma while our computational method 
is compatible with it. It shall be noted, however, that 
the recent observation of broken parity and time reversal 
evident in polarized optical images of planar chiral struc- 
tures Q| calls for re-examining the validity of the Lorentz 
lemma for planar chiral structure. 
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